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Patterning of dielectric nanoparticles using dielectrophoretic forces generated 
by ferroelectric polydomain films 
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Czech Republic 
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A theoretical study of a dielectrophoretic force, i.e. the force acting on an electrically neutral particle in 
the inhomogeneous electric field, which is produced by a ferroelectric domain pattern, is presented. It has 
been shown by several researchers that artificially prepared domain patterns with given geometry in fer- 
roelectric single crystals represent an easy and flexible method for patterning dielectric nanoobjects using 
dielectrophoretic forces. The source of the dielectrophoretic force is a strong and highly inhomogeneous 
(stray) electric field, which exists in the vicinity of the ferroelectric domain walls at the surface of the ferro- 
electric film. We analyzed dielectrophoretic forces in the model of a ferroelectric film of a given thickness with 
a lamellar 180° domain pattern. The analytical formula for the spatial distribution of the stray field in the 
ionic liquid above the top surface of the film is calculated including the effect of free charge screening. The 
spatial distribution of the dielectrophoretic force produced by the domain pattern is presented. The numerical 
simulations indicate that the intersection of the ferroelectric domain wall and the surface of the ferroelectric 
film represents a trap for dielectric nanoparticles in the case of so called positive dielectrophoresis. The effects 
of electrical neutrality of dielectric nanoparticles, free charge screening due to the ionic nature of the liquid, 
domain pattern geometry, and the Brownian motion on the mechanism of nanoparticle deposition and the 
stability of the deposited pattern are discussed. 



I. INTRODUCTION 

A well-known fact follows from the Maxwell theory 
that the inhomogeneous electromagnetic field induces a 
force even on electrically neutral particles. However, for 
a long time there has been no attention on how this phe- 
nomenon could be practically applied. A new impetus in 
this field was provided in a 1951 paper by Pohl^ where he 
studied the influence of inhomogeneous electrostatic field 
on the motion of small plastic particles suspended in in- 
sulating dielectric fluids. Pohl named this phenomenon 
dielectrophoresis (DEP) and it was shown in his 1978 
monograph^ that DEP forces can be usefully used for 
the manipulation of small particles. This book, in fact, 
motivated studies of other possibilities to use the elec- 
tric field for trapping or controlling the motion of small 
particles. 

At the same time, in 1980s, a new field of nan- 
otechnology emerged and DEP has become an impor- 
tant and intensively studied tool for the manipulation 
with nanoparticles. Examples to be mentioned here 
are the DEP assembly of nanowires from nanoparti- 
cle suspensions^, DEP precise positioning of carbon 
nanotubes^, particle separation using DEP^"^, DEP for- 
mation of nanodroplets^ or nanofibcrs^, etc. In all the 
aforementioned applications, the inhomogeneous electric 
field, which is essential for the existence of DEP forces, 
is achieved by properly arranged electrodes. The cru- 
cial point related to nanotechnology is that, when there 
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is a wish to work with nanoparticles, it is necessary to 
deposit very thin electrodes which are very close to each 
other. This task can be obviously managed^" but, in fact, 
it is usually technologically demanding. In addition, a 
fixed set of electrodes unfortunately does not allow flexi- 
ble changes in the nanoparticle pattern geometry, which 
may be necessary in some advanced applications. 

In order to overcome disadvantages of the fixed elec- 
trode approach, attention of researchers has been drawn 
to answer the question: would it be more useful to use 
an inhomogeneous field existing above the surface of a 
ferroelectric film with a domain pattern? The interest 
of scientists in this question has been stimulated by a 
great advance that has recently been achieved in a pre- 
cise writing of ferroelectric domains of nanometer dimen- 
sions. Several methods for nanoscale domain control in 
ferroelectric films have been developed by several groups 
of researchers. Classical examples of domain engineer- 
ing to be mentioned here are the domain poling with 
patterned electrodes^"'^^^'^, domain patterning by scan- 
ning probe microscopy and focused electron or ion 
beam domain engineering^^"^^. Details about advanced 
nanoscale domain engineering methods can be found in 
the recent review paper by Li and BonnelP^. The afore- 
mentioned methods make it possible to prepare arbitrary 
ferroelectric domain patterns in a flexible way. 

In fact, there exist two different means to use ferro- 
electric domains for nanoparticles patterning. The first 
approach is based on the use of domain-specific chemi- 
cal reactions and it is called ferroelectric lithography^^ 
Using this method, the direct assembly of virus particles 
has been realized by Dunn et a/.^^The second approach 
is based on the use of physical DEP forces that are pro- 



duced by bound charges at the surface of the ferroelectric 
film with a domain pattern, which has been recently re- 
ported by Ke^^ et al. and Grilli et al?'' . 

It is the second approach that has motivated the the- 
oretical analysis presented below, where we will address 
the study of the DEP forces produced by lamellar 180° 
domain pattern. In Sec. II, we first present the details 
of our model of a ferroelectric film with 180° domain 
domain pattern. Section III presents a straightforward 
calculation of the electrostatic potential and the electric 
field produced by the ferroelectric domain pattern. We 
will cover the general configuration of the domain pat- 
tern including the effect of free charge screening as well 
as two cases of a special interest. First, the domain pat- 
tern with wide domains compared to the thickness of the 
ferroelectric film (Sec. Ill A) and the system with a dense 
domain pattern (Sec. IIIB). In Sec. IV, calculation and 
numerical simulations of the DEP forces produced do- 
main patterns in four systems with different configura- 
tions are presented and compared. Section V presents 
a brief discussion of the effect of electrical neutrality of 
dielectric nanoparticles on the mechanism of nanopat- 
terning (Sec. VA), the effect of free charge screening 
on the inhomogeneous field and the value of DEP force 
(Sec. VB), the role of domain pattern geometry on the 
value of the DEP force produced by the ferroelectric do- 
main pattern (Sec. VC), the role of the Brownian motion 
on the stability of the nanoparticle patterns (Sec. VD), 
and, finally, the estimation of the values of important nu- 
merical parameters that are essential in the design of a 
ferroelectric system used for patterning of nanoparticles 
(Sec. VE). 



II. GEOMETRY AND MATERIAL PROPERTIES OF 
THE FERROELECTRIC FILM 

Figure 1 shows the geometry of the considered model 
of a ferroelectric film with a lamellar 180° ferroelectric 
domain pattern crossing the film thickness d. The space 
above the top surface of the film is filled with an ionic 
liquid of permittivity Em and with a given concentration 
of free charge carriers. The bottom surface of the fer- 
roelectric film is attached to the ground electrode. We 
consider that the top and bottom surfaces of the film 
are perpendicular to the ferroelectric axis z of the at- 
tached Cartesian coordinate system. The axis x of the 
coordinate system is perpendicular to the system of 180° 
domain walls. Using the "hard ferroelectric" approxima- 
tion, we express the electric displacement of the ferro- 
electric as a sum of the linear dielectric response of the 
crystal lattice to the electric field and the constant spon- 
taneous polarization Pq (whose orientation differs from 
domain to domain): 
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FIG. 1. Geometry of the ferroelectric film of the thickness 
d with the 180° domain pattern. Gray arrows indicate the 
orientation of the vector of spontaneous polarization in the 
film. Dashed lines indicate positions of the 180° domain walls. 
Symbols a+ and a_ stand for the widths of domains where the 
vector of the spontaneous polarization is oriented along and 
against the orientation of the z-axis of the attached coordinate 
system, respectively. Symbols Sa and Ec stand for the relative 
permittivity tensor components of the ferroelectric, symbol 
Em stands for the permittivity of the liquid. 



where £c and Sa are the components of the permittivity 
tensor of the crystal lattice in the directions parallel and 
perpendicular to the ferroelectric axis, respectively, and 
£o is the vacuum permittivity. We consider that the vec- 
tor of spontaneous polarization is perpendicular to the 
film surfaces and its magnitude Pq is constant within 
each domain. For given widths a+ and a_ of domains, 
where the vector of the spontaneous polarization is ori- 
ented along and against the direction of the z-axis of the 
attached coordinate system, respectively, we define the 
average (net) spontaneous polarization of the ferroelec- 
tric layer P/v = Po (a+ — a_)/(2a), where a = (a++a_)/2 
is the domain spacing. 

The abrupt change of the spontaneous polarization to 
zero at the interface of the ferroelectric film and the ionic 
liquid yields the appearance of a bound charge of the sur- 
face density ah{x), which is the source of the spatially 
nonuniform electric field. The spatial distribution of the 
bound charge due to the discontinuous change of the nor- 
mal component of the spontaneous polarization at the 
ferroelectric film surface in the direction of the x axis 
can be written in the form: 



ab{x)^PN+} 

' ^ ■wn 



Dx 



SoScEz ± Po, 



(la) 
(lb) 



where k — ix j a. 



sm 



n-K ( Pm 
2 I Po 



cosnkx, (2) 



3 



III. ELECTROSTATIC POTENTIAL PRODUCED BY 
THE DOMAIN PATTERN 

In this section we present the calculation of the elec- 
trostatic potential produced by the ferroelectric domain 
pattern. Since the dielectric properties of the ferroelec- 
tric film and the ionic liquid are quite different, we use 
the following symbols and ipj for the electrostatic po- 
tential in the ionic liquid (for z > 0) and the ferroelectric 
film (for > > —d), respectively. In our analysis, 
we consider the general situation where the ionic liquid 
above the top surface of the ferroelectric film contains 
free charge carriers of given concentration. In such a sit- 
uation, the electric field is screened by the free charge car- 
riers and vanishes within some typical distance from the 
surface of the ferroelectric film. In order to calculate the 
spatial distribution of the electrostatic potential in the 
system, we adopt the Debye-Hiickel theory^® (DHT) and 
the aforementioned functions for the electrostatic poten- 
tial should satisfy: first, the following partial differential 
equations of DHT and electrostatics 
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where 



2NAe^I 



(3a) 
(3b) 

(3c) 



is the Debye screening length, ks is the Boltzmann 
constant, T is the thermodynamic temperature, Na is 
the Avogadro number, e is the electron charge, / ~ 
ionic strength of the liquid, where 
Ci and Zi are the molar concentration and the charge 
number of the i-th particular ion in the liquid. 

Second, the boundary conditions for the continuity of 
the electrostatic potential 
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iff 
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Third, the boundary condition for the continuity of the 
tangential component of the electric field 
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at z = 0, 
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Fourth, the boundary condition for the continuity of the 
normal component of the electric displacement 
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' dz 
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It is a straightforward task to show that the the follow- 
ing functions in the form of Fourier series represent the 



solution of the electrostatic problem given by Eqs. (3): 
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where x = kx, z = kz, p = ttPn / {2Po), c ~ \J Eajsc-, 
g = cEc/Sm, R = ckd, and 
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In the special case of the absence of free charge carriers 
in the liquid, i.e. infinite value of the Debye screening 
length Xu, the values of parameters and rjuin) for 
all n > 1 tend to 1 and the solution of the electrostatic 
problem given by Eqs. (3) reduces down to the form that 
is already available in literature^^ , however in the special 
case of the ferroelectric film in a vacuum, i.e. for = 1. 

Equations (4) represent the exact solution of the elec- 
trostatic potential produced by the domain pattern. For 
the sake of clearer graphical presentation, it is convenient 
to use the normalized value $ of the electrostatic poten- 
tial: $ = tp{k/CE)- In the following subsections, we will 
present two special situations in the system with the ab- 
sence of free charge carriers, i.e. Au 3> a, in which the 
above expressions can be further simplified. 



A. Systems with thin ferroelectric films and wide domains 

The first situation is considered for the systems with 
an isolating dielectric liquid where the thickness of the 
ferroelectric film d is smaller than the average distance 
between the domain walls a, i.e. i? <C 1. In addition, it 
is reasonable to consider that the permittivity of ferro- 
electric film is much larger than the permittivity of the 
dielectric liquid, i.e. g ^ 1. In this case, the factors 
(1-1- g coth nR) in the denominator of the leading terms 
in Eq. (4a) are dominated by g coth nR and it is quite ac- 
ceptable to approximate {1+ g coth nR) by ^(l-l-cothni?). 
Using the Euler formula e*^ = cosz 4- isinz, using the 
sum and difference formulas for the trigonometric func- 
tions, and using the following formula: 
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FIG. 2. Contour plot of the elec- 
trostatic potential produced by the 
domain pattern in the system with a 
thin ferroelectric film and wide do- 
mains, 7? = 0.1, (a, c) and with 
a dense domain pattern, R — 10, 
(b, d). In addition, situations with 
the absence of free charge carriers, 
kXo = 10^, (a, b), and with free 
charge screening, kXo = 0.2, (c, 
d), are compared. The electrostatic 
potential is normalized to the con- 
stant (CE/k) given by Eq. (4c). Red 
arrows represent the inhomogeneous 
electric field produced by the do- 
main pattern. The following values 
of material parameters were consid- 
ered p = 0.5, g = 15.8. Note the dif- 
ference in the peak values of the nor- 
malized electrostatic potential and in 
the normalized electrostatic poten- 
tial gradient in the vicinity of the do- 
main walls in the cases (a, c) and (b, 
d). 



d) ii = 10 and kXo = 0.2 
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where Lig z = z^/n'^ is the polylogarithm function, 

the electrostatic potential in the dielectric medium can 
be expressed in a form: 



k I a ' a ' 2a V £c ' 2Po ' ' ^ ' 



where the function stands for the sum in Eq. (4a) 
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with wide domains. It should be noted that in the case 
of systems where the permittivity of the dielectric liq- 
uid, Em, is approximately equal to y/ec£a, i-e. g « 1, 
the above expressions Eqs. (6) equal exactly those in 
Eq. (4a). 

The electric field in the liquid produced by the ferro- 
electric domain pattern with wide domains can be ex- 
pressed in a straightforward way: 
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The superscript '^"'^ indicates that the functions are ex- 
pressed in the approximation of the thin ferroelectric film 



dz 



N 



r s"^*") ( — — — ^-Pj 

' 1 ' a ' 2a V Sc ' 2Po 



(7b) 



where 



<^'f\x,Z,R,p) - -^lo; 



1 [cosh z — sin(a; — p)] [cosh(z -I- 2R) + sm{x + p)] 



[cosh z + sin(i; + p)] [cosh(z -f 2R) — sm{x — p)] ' 



(8a) 



5 



(,,,\ , ^ 1 2 cos i5 (e^ cos £ + sin io) ^ ^ 

^rUx,z,R,p) = -— arctan , , , „ . + 8b 

4g (e^ + 2 cosx smp) — cos2p 

1 2cosp (cosi + e~^~^^ sinp) 
■ arctan ■ 



Ag 2 sinp [cos x + cosh(z + 2R) sinp] + 2 cos^ p sinli(z + 2R) 



In the derivation of Eqs. (8a), we have used the property 
of the polylogarithm function d{lj\2 z) / dz ~ — log(l — 
z)/z and the definition of the logarithm in a complex 
plane log z = log | ^ | + i arg z . 



B. Systems with dense domain patterns 

The second situation, for which it is possible to write 
down a closed formula for the electrostatic potential and 
the components of the electric field, is the ferroelectric 
film with a dense domain pattern, i.e. the systems where 
the thickness of the ferroelectric film d is larger than the 
average distance between the domain walls a, i.e. i? ^ 1. 
In this case, the terms coth nR in the denominator of the 
leading terms in Eq. (4a) can be replaced by 1 and, thus, 
it is quite acceptable to approximate (1 + gcothni?) by 
(1 +5). Following the same algebraic procedures as in 
the previous subsection, the electrostatic potential in the 
dielectric liquid can be expressed in a form: 



k V a ' a ' 2Po 
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where the function ^^'^^ stands for the sum in Eq. (4a) 
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The superscript ('^^ indicates that the functions are ex- 
pressed in the approximation of the dense domain pat- 
tern. 

The electric field in the dielectric liquid produced by 
the dense ferroelectric domain pattern can be expressed 
in a straightforward way: 
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C. Simulation of the electrostatic potential and 
electrophoretic force 

Figures 2a and 2b show the spatial distribution of the 
electrostatic potential (p given by Eq. (4a) and the vec- 
tors of the electric field produced by the periodic domain 
pattern in the dielectric liquid in the two aforementioned 
situations characterized by the absence of the free charge 
carriers in the liquid medium (i.e. with a large value of 
the Debye screening length kXn = lO'^). Figure 2a shows 
the system with a thin ferroelectric film and wide do- 
mains, R = 0.1, and Figure 2b shows the system with 
a dense domain pattern, R = 10, (b). The electrostatic 
potential ip is normalized to the constant {CE/k) given 
by Eq. (4c) and, thus, the contour plots correspond to 
the spatial distributions of the dimensionless functions $ 
given by Eq. (6b) and (9b). Figures 2c and 2d present 
the effect of the inhomogeneous electric field screening 
by free charge carriers in the ionic liquid (i.e. the small 
value of the Debye screening length compared to the do- 
main spacing, kXjj = 0.2) in the systems with the same 
geometry as in Figs. 2a and 2b. In the bottom part of 
the figure, there are indicated the positions of the 180° 
domain walls and the thickness of the ferroelectric film. 
Red arrows represents the vectors of the inhomogeneous 
electric field produced by the domain pattern and they 
are proportional to the electrophoretic (EP) force, i.e. 
the Coulomb force acting on charged particles. 

If we consider an approximately spherical particle of 
a diameter r with free charges on its surface of a sur- 
face density cr, the electrophoretic force acting on such a 
particle is given by the classical formula: 



(12) 



With use of the notation adopted in this Article, Eq. (12) 
can be rewritten in the form: 



where 



(13a) 



(13b) 



and the symbol V stands for the gradient operator in the 
ii-coordinate system, where Xi = kxi. 

In the numerical calculation of the normalized electro- 
static potential, we considered the values of the ionic liq- 
uid permittivity Era = 3 and the in-plane and out-of plane 
permittivity of the ferroelectric Ea = 84 and Ea = 29, re- 
spectively, which corresponds to the dielectric properties 
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FIG. 3. Contour plot of the pos- 
itive DEP force potential produced 
by the domain pattern in the sys- 
tem with a thin ferroelectric film and 
wide domains, R = 0.1, (a, c) and 
with a dense domain pattern, R — 
10, (b, d). In addition, situations 
with the absence of free charge car- 
riers, kXo = 10^, (a, b) , and with 
free charge screening, kXo = 0.2, 
(c, d), are compared. The potential 
is normalized to Cdep- Red arrows 
represents the DEP force produced 
by the domain pattern. The fol- 
lowing values of material parameters 
were considered p = 0.5, g = 15.8. 
Note the difference in the peak values 
of the normalized dielectrophoretic 
force potential in the vicinity of the 
domain walls in the cases (a, c) and 
(b, d), which is controlled only by the 
geometry of the domain pattern. 



0.2 



of lithium niobate'^^''^^. The three aforementioned values 
of permittivity give the value of the parameter g — 15.8. 
Finally, the non-neutral domain pattern was considered 
with p = 0.5, which corresponds to Pn/Pq = 0.32. It 
is noticeable to see the difference in the peak values of 
the normaUzed electrostatic potential -0.005 and 0.005 in 
Fig. 2a; and 0.197 and 0.289 in Fig. 2b, which is not con- 
trolled by the values of spontaneous polarization Pq but 
by the ratio of the ferroelectric film thickness over the 
domain spacing d/a, which is proportional to the value 
of the parameter R. The similar effect of the value of 
R can bee seen on the normalized electrostatic potential 
gradient in the vicinity of the domain walls in Figs. 2a 
and 2b. 



IV. DIELECTROPHORETIC FORCE 

The dielectrophoretic force Fdep that acts on approxi- 
mately spherical particle of a diameter r and the dielec- 
tric constant £p in the inhomogeneous electric field E is 
given by the classical formula^: 

Fdep = KdepV {El + El) , (14a) 

where 

Kdep = 27rr3£,„£o ■ (Mb) 

Sp -r ZiSjyi 

It should be noted that depending on the permittivity 
values of the medium and particles, the situation char- 
acterized by Ep > £m, i.e. Kdep > 0, is called positive 



dielectrophoresis and the situation, when Sp < i.e. 
Kdep < 0, is called negative dielectrophoresis. 

With use of the notation adopted in this Article, 
Eq. (14a) can be rewritten in the form: 

Fdep = KdepClW{<i>l+<^j) (15) 
= Kdep elk V {<i>l+^j) . 

For the purpose of the graphical presentation of spatial 
distribution of the DEP forces generated by the domain 
pattern, it is convenient to introduce the DEP force po- 
tential by formula Fdep = ^VC/^ep- Thus, the function 
Udep can be expressed in the form 

Udep = -Cdep ($1 + $1) , (16a) 

where 

— 32 Pq (gp - £m) 
Cdep - — — — , „ ^ ■ (16b) 

a to £m \£p + ^fcmj 

Figures 3a-d show the spatial distribution of the di- 
electrophoretic field potential Udep given by Eq. (16a) 
and the vectors of the dielectrophoretic force produced 
by the periodic domain pattern in the ionic liquid in the 
same four situations as in Figs. 2. Figures 3a and 3c 
show the systems with a thin ferroelectric film and wide 
domains, R = 0.1. On the other hand. Figures 3b and 3d 
show the systems with a dense domain pattern, R = 10. 
The dielectrophoretic field potential Udep is normalized to 
the constant Cdep given by Eq. (16b) and, thus, the con- 
tour plots correspond to the spatial distributions of the 
dimensionlcss functions — ($| -|- <&?). Figures 3a and 3b 
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show the systems with a dielectric hquid and a low con- 
centration of free charge carriers, which corresponds to 
the large value of the Debye screening length, kXo = 10"^. 
On the contrary. Figs. 3c and 3d show the systems, where 
the inhomogeneous electric field is screened due to free 
charge carriers in the ionic liquid and where the Debye 
screening length is smaller than the domain spacing, i.e. 
kXo = 0.2. It is noticeable to sec the difference in the 
peak values of the normalized dielectrophoretic field po- 
tential in Figs. 3a and 3c, and 3b and 3d, respectively. 
Again the point is that the peak values are not controlled 
by the values of spontaneous polarization Pg but by the 
ratio of the ferroelectric film thickness over the domain 
spacing d/a, which is proportional to the value of the 
parameter R. 



V. DISCUSSION AND CONCLUSIONS 

Theoretical results presented in the above Sections in- 
dicate qualitative features that play the key role in the 
decoration of nanoparticles using polydomain ferroelec- 
tric films. In the following subsections, we will briefly 
discuss (i) the effect of electrical neutrality of dielectric 
nanoparticles on the mechanism of nanopatterning, (ii) 
the effect of free charge screening on the inhomogeneous 
field and the value of DEP force, (iii) the role of domain 
pattern geometry on the value of the DEP force pro- 
duced by the ferroelectric domain pattern, (iv) the role 
of the Brownian motion on the stability of the nanoparti- 
cle patterns, and, finally, (v) we will estimate the values 
of important numerical parameters that are essential in 
the design of a ferroelectric system used for patterning 
of nanoparticles. 



A. The role of electroneutrality of dielectric nanoparticles 

Since the electrical neutrality of nanoparticles is rarely 
the case in real experimental situations, it is necessary 
to compare the EP and DEP forces and to analyze the 
conditions that controls switching between the EP and 
DEP mechanisms of nanoparticle deposition. 

Figure 4 shows the comparison of the EP and DEP 
force dependences on the ratio x/a'm the system with a 
thin ferroelectric film and wide domains, i? = 0.1 (solid 
line), and in the system with a dense domain pattern, 
i? = 10 (dashed line), at the distance z = 0.01 a from the 
surface of the ferroelectric film. The EP force value is 
normalized to C^p given by Eq. (13b). DEP force value 
is normalized to Cdep given by Eq. (16b) and, therefore, 
corresponds to the positive DEP force. It is seen that in 
the case of the EP mechanism, the x-component of the 
EP force is zero in the middle of each domain and the 
potentially charged nanoparticles would be preferentially 
deposited in the middle of those domains that have the 
opposite charge than the charged nanoparticle. On the 
contrary, the x-componcnt of the DEP force is zero at 
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FIG. 4. Comparison of the a::-component of the EP and DEP 
force dependences on the ratio a;/a in the system with a thin 
ferroelectric film and wide domains, R = G.l (solid line), and 
in the system with a dense domain pattern, 7? = 10 (dashed 
line), at the distance z = 0.01 a from the surface of the ferro- 
electric film. The EP force value is normalized to Cep given 
by Eq. (13b). DEP force value is normalized to Cdep given 
by Eq. (16b) and, therefore, corresponds to the positive DEP 
force. Note the qualitative difference between the EP and 
DEP force spatial distribution. 



the domain wall position, abruptly increases with the in- 
creasing distance from the domain wall and reaches its 
maximum F^^max at the distance \xmax from the do- 
main wall. It is seen that the reduction of the thickness 
of the film by a factor of 100, the peak values of the DEP 
force x-component a smaller by a factor of 10. It is essen- 
tial that the area of the intersection of the domain wall 
with the surface of the ferroelectric film represents a trap 
for nanoparticles in the case of a positive dielectrophore- 
sis. 

Figure 4 shows that the peak normalized values of the 
EP and DEP forces are of the same order of magnitude. 
Therefore, one should compare the parameters Cep and 
Cdep to get the critical value of the surface charge that 
controls the switching between the EP and DEP nanopar- 
ticle deposition mechanism. As a result, the DEP mech- 
anism of deposition is dominant over the EP mechanism, 
when the surface charge density of the nanoparticles is 
much smaller than the critical value: 



2Por (gp - £m) 
a (2e,n + e„) 



(17) 



Table I presents numerical values of materials param- 
eters considered in our theoretical analysis for lithium 
tantalate (LiTaOs), lithium niobate (LiNbOs), barium 
titanatc (BaTiOs), and lead titanate (PbTiOa) that can 
be prepared in a form of single crystal samples. In our 
numerical estimates, we further adopt parameters of ma- 
terials that have been used previously in experiments 
by Ke et al.^^, In particular, we consider the use of 
lithium niobate single crystal with the domain spacing 
a = 5 /im, polystyrene nanoparticles with dielectric con- 
stant Sp — 2.7 and radius r = 65 nm, and dodccane liquid 
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TABLE I. Material parameters (spontaneous polarization Po, 
permittivity tensor components in the direction parallel and 
perpendicular to the direction of the spontaneous polarization 
Ec and Sa, respectively, and the pyroelectric coefficient 7s 
at the room temperature) of selected ferroelectric materials: 



LiTa03^^'^^ LiNbOa^"'^^ 



BaTiOa 



and PbTiOa^^ 



Parameter 


Unit 


LiTaOs 


LiNbOa 


BaTiOg 


PbTiOs 


Po 


Cm-2 


0.55 


0.75 


0.26 


0.76 




1 


44 


29 


188 


66 


£a 


1 


53 


84 


3600 


124 


7s 


/iCm-^K"^ 


230 


40 


200 


380 



with £,„ = 2, the critical value of the nanoparticlc surface 
charge is acrit ~ 0.002 C m"^. 



B. The role of free charge screening 

Since the nanoparticlc deposition experiments are car- 
ried out in liquid media, the presence of free charge carri- 
ers in the liquid cannot be completely avoided and some- 
times it is even induced intentionally. It means that the 
effect of finite Debye screening length on the values of EP 
and DEP forces should be carefully analyzed. At first let 
us consider the simplest case of water. At room temper- 
ature, the permittivity of water is about 80 and the ionic 
strength equals / = 10^** molm~'^, which yields the typi- 
cal value of the Debye screening length A^j « 1 /.im. In the 
case of nonpolar dielectric liquids, such as pentanoic acid 
or silicone oil, the value of the Debye screening length 
can be of several orders of magnitude larger. 

Figure 2 clearly shows that the inhomogeneous elec- 
tric field essentially decays within the distance of about 
2a from the surface of the ferroelectric film. This means 
that the decoration of nanoparticles using DEP forces 
should not be affected by free charge carriers in the sys- 
tems where the Debye screening length is larger than the 
domain period 2a. From this condition, one can calculate 
the critical ionic strength of the lic^uid medium: 



KIc 



£o£mkBT 



(18) 



Considering the experiments by Grilli^^ with pentanoic 
acid Em = 2.66 and the value of domain spacing 
a = 100 //m, the Eq. (18) yields the value Icrit = 
T.riO^'^molm-^. 



C. The role of domain pattern geometry 

Both Figs. 2 and 3 clearly indicate an essential role 
of the geometry of the domain pattern in the numeri- 
cal value of the dielectrophoretic force produced by the 
bound charges of the surface of the ferroelectric film with 
a domain pattern. To get a clear insight into the role of 



domain pattern geometry on the DEP force. Fig. 5 shows 
the maximum value of the x-component of the DEP force 
and the position of this maximum, which is expressed 
in the distance from the domain wall and measured as 
a fraction of the domain spacing a, were calculated as a 
function of the parameter R for three distances z from the 
surface of the ferroelectric film. Figure 5 indicates that 
there is a qualitative difference in the aforementioned de- 
pendencies for values of R smaller and greater than one. 
In the systems with a thin ferroelectric film and wide do- 
mains, which are characterized by i? < 1, the maximum 
value of the a;-component of the DEP force is decreasing 
polynomially with a decreasing value of the parameter R. 
On the other hand, the systems with a dense ferroelec- 
tric domain pattern, which is characterized by i? > 1 , the 
maximum of the DEP force is insensitive to the value or 
R. It should be noted that the maximum values of the 
DEP force x-component at the given distance z = z/k 
from the surface of the ferroelectric film are normalized 
to Cdep and they are, therefore, dominantly controlled by 
the geometry of the domain pattern. 



D. The role of Brownian motion 

In order to achieve a stable pattern of nanoparticles, 
the DEP force that acts to move the nanoparticlc into 
the trap must overcome the forces that acts on that par- 
ticle due to thermal movements of surrounding medium 
molecules, i.e. the forces that are responsible for the 
Brownian motion. The role of Brownian motion in DEP 
trapping of particles has been already analyzed by sev- 
eral authors'^^''^''. Here we adopt the results achieved by 
Hughes and Morgan'^^ where the limit for particle trap- 
ping is estimated by comparing the velocities that are 
induced by DEP force with that of the Brownian mo- 
tion. Result of this analysis yield the formula for the 
estimated critical minimal particle radius that is needed 
for successful trapping of a particle:'^^ 



r > Trrit = 



lOkeT 



Trera£oKAL\VE^\ 



1/3 



where 



2er, 



(19) 



(20) 



is the Clausius-Mossotti factor, AL is typical distance 
over which the field gradient Vi?^ has a given minimal 
value (see the inset in Fig. 5b), fc^ is the Boltzmann 
constant and T is the thermodynamic temperature. In 
the notation of this Article, the above formula can be 
rewritten in the form: 



lOfcsT 



■nem£QKC'^S{z,g,R) 



1/3 



where 



5(z,g,i?) = AL ($l + 1>i) 



(21) 



(22) 
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FIG. 5. Dependence of the s-component of the positive DEP force maximum value normalized to Cdep (a) and the distance 
Axmax of the point, where the x-component of the DEP force reaches its maximum value (b), from the domain wall (normalized 
to the domain spacing o) as a function of the parameter R = {-nd) / {2a)^ea/ £c- The inset shows the maximum value of the 
DEP force Fx,max/Cdep at the distance l^Xmax from the domain wall and the distance AL over which the DEP force has a 
given minimal value. The parameters /^Xmax and AL are measured as a fraction of domain spacing a. 



is a dimcnsionless function given by the product of the ef- 
fective value of the function | Vi(<i>|-|-$|)| on the distance 
AL normahzed to the domain spacing a, i.e. AL = k AL. 
For a given material, given values of the parameters g 
and i?, and a given distance from the surface of the fer- 
roelectric film z, it is possible to find a maximum value 
of the function S'max = 'niax.{S {z , g , R)} by a numerical 
calculation (see Tab. II). 



E. Conclusions 

The spatial distribution of the inhomogencous electric 
field produced in a liquid above the top surface of the 
poly domain ferroelectric film has been calculated. The 
inhomogeneous electric field is the source of the elec- 
trophoretic and dielectrophoretic forces that can be used 
for a deposition of nanoparticle pattern. The key element 
of our analysis, considering the ionic nature of the liquid, 
makes our model quite realistic. Our numerical simula- 
tions have shown that the intersection of the domain wall 
with the surface of the ferroelectric film represents a trap 
for dielectric nanoparticles in the case of positive dielec- 
trophoresis. The analysis of the critical conditions that 
control the deposition of dielectric nanoparticles and the 
stability of the deposited pattern are presented. 

As an application of our analysis, one can consider two 
ways of realization of DEP force patterning using ferro- 
electric domains. First, it is the situation where the in- 
homogeneous electric field is produced by a spontaneous 
polarization reversal to form the required domain pat- 
tern using advanced methods of the nanoscalc domain 
geometry engineering. Considering the values of mate- 
rial parameters that correspond to experiments by Grilli 
et al.^"^ , i.e. pentanoic acid with e,„ = 2.66 and flour 
particles with Sp = 3, the estimated numerical parame- 
ters of positive DEP force are presented in Tab. II. It is 



TABLE II. Calculated characteristic numerical parameters of 
the DEP force produced due to spontaneous polarization re- 
versal from a single domain state to a polydomain pattern 
of a given geometry in the selected ferroelectric materials. It 
is considered: r = 1 /xm, z — 0.01 a, R — 10, Em ~ 2.66, 
£p = 3.0, and T = 300 K. 



Parameter 


Unit 


LiTaOs 


LiNbOa 


BaTiOs 


PbTiO 


9 


1 


18.1 


18.6 


309.3 


34.0 


Ce 


GVm-^ 


30 


41 


14 


41 


Cdep 




168 


312 


37.5 


321 


1 V-E''^ |max 


lO^^V^m-^ 


110 


195 


0.09 


60 


-P^dep, max/ Cdep 


10^^ 


66 


63 


0.25 


20 






8.1 


14.4 


0.007 


4.6 


5'max 


10"^ 


5.39 


5.14 


0.02 


1.6 


l^crit 


nm 


1.5 


1.2 


15.8 


1.8 



seen that using this procedure, it is possible to achieve 
a large local values of the DEP forces that results in 
very small critical radii of particles (about 2 nm) that 
can be trapped near the domain walls. Disadvantage of 
this method stems from the fact that the bound charges 
due to the discontinuous change of the spontaneous po- 
larization at the surface of the ferroelectric film can be 
partially or completely compensated by free charge car- 
riers due to the possible nonzero conductivity of the fer- 
roelectric film or due to the ionic nature of the liquid. In 
this case, the inhomogeneous electrostatic field can dis- 
appear after some time and other method nanoparticle 
decoration should be used. 

The second way of the realization of DEP forces using 
ferroelectric polydomain sample is the use of pyroelectric 
effect. It is known that the pyroelectric coefficients in the 
opposite domains are of the opposite sign. Therefore, 
when the ferroelectric polydomain sample is heated up 
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TABLE III. Calculated characteristic numerical parameters 
of tire DEP force produced due to pyroelectric effect after 
heating the ferroelectric film from the room temperature up 
by AT = 50K in the selected ferroelectric materials. 



Parameter 


Unit 


LiTaOa 


LiNbOs 


BaTiOa 


PbTiOs 


Ce 


GVm-i 


0.62 


0.11 


0.54 


1.1 




nN 


73.4 


2.2 


55.5 


222 


1 V-E"^ max 


lO^'^V^m-^ 


48 


1 


0.1 


40 




nN 


3.6 


0.1 


0.01 


3.2 




nm 


20 


64 


140 


20 



by AT, the bound charges of the surface density ±75 AT 
appear on the surface of the film due to the pyroelec- 
tric effect, where 75 is the pyroelectric coefficient and 
its sign differs from domain to domain. In order to get 
some numerical estimates for this case, the numerical val- 
ues of parameters Ce and Cdep should be calculated by 
replacing the symbol Pq by 7sAT in Eqs. 4c and 16b. 
Estimated numerical parameters of DEP forces achieved 
by this method are presented in Tab. III. Advantage of 
this approach is its relatively easy way of realization, e.g. 
using local laser beam heating of a single domain film or 
uniform heating of a polydomain film. Disadvantage is 
a presence of smaller values of DEP forces that result in 
greater value of the critical radius of particles that can be 
trapped at required locations (larger than 20 nm). The 
numerical estimates presented in Tab. Ill are in a qual- 
itative agreement with experimental results by Grilli et 

Finally, we have theoretically demonstrated that using 
ferroelectric domain patterns it is possible to achieve a 
very flexible way of nanoparticle decoration, which has 
several advantages: (i) using advanced and recently well 
developed methods of domain geometry engineering, it is 
possible to achieve a "rewritable" systems for the prepa- 
ration of an arbitrary nanoparticle patterns, (ii) since the 
ferroelectric domain wall is an extremely thin spacial fea- 
ture comparing with other ferroic material or deposited 
electrodes, it is possible to realize patterning of nanopar- 
ticles with sub-micron lateral dimensions. We believe 
that the obtained results represent a useful tool that can 
be profitably used for designing the systems with ferro- 
electric polydomain films representing a flexible way for 
decoration of nanoparticles. 
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